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TRAVELING WAVES FOR THE MASS IN MASS MODEL OF 

GRANULAR CHAINS 

PANAYOTIS G. KEVREKIDIS, ATANAS G. STEFANOV, AND HAITAO XU 


Abstract. In the present work, we consider the mass in mass (or mass with mass) 
system of granular chains, namely a granular chain involving additionally an internal 
resonator. For these chains, we rigorously establish that under suitable “anti-resonance” 
conditions connecting the mass of the resonator and the speed of the wave, bell-shaped 
traveling wave solutions continue to exist in the system, in a way reminiscent of the results 
proven for the standard granular chain of elastic Hertzian contacts. We also numerically 
touch upon settings where the conditions do not hold, illustrating, in line also with recent 
experimental work, that non-monotonic waves bearing non-vanishing tails may exist in 
the latter case. 


1. Introduction 

A topic that has been progressively gaining more attention over the last two decades 
within the nonlinear dynamical systems of the famous Fermi-Pasta-Ulam type [9], concerns 
the study of the so-called granular crystals [25l|28l[l5l|32]. The latter consist of chains of 
elastically interacting (through the so-called Hertzian contacts) beads that are not only very 
experimentally accessible, but also extensively tunable and controllable, as regards their 
materials, geometry, heterogeneity, etc. Another reason for the considerable appeal of such 
simple lattice structure, which can be tailored to be one-, two- or even three-dimensional, 
is the wealth of nonlinear excitations that have arisen in such settings, and which include 
robust traveling waves, bright and dark breather structures, and shock waves among others, 
as summarized in the above reviews. Finally, yet another element promoting the interest 
in such setups concerns their potential relevance to a broad range of applications such as 
actuating devices nzi, acoustic lenses [29], mechanical diodes [221I2I1I2], logic gates [20] 
and sound scramblers [3 [2S] • 

Arguably, the most significant and well-studied excitation, the principal workhorse 
around which many theoretical analyses and experimental results have been centered in 
such granular systems is the traveling wave. This structure was originally established in the 
pioneering work of Nesterenko in [27] (see also the review of [25]). While it was originally 
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believed to be genuinely compact, asymptotic analysis illustrated its doubly exponential 
nature (in the granular system without precompression; in the presence of precompres¬ 
sion the decay of the wave becomes exponential) [HI [1]. /^From a mathematically rigorous 
perspective, the work of [lO] enabled as a special case example application the proof of ex¬ 
istence of such a traveling wave in the work of [ 23 ] , while later the framework of [8] (earlier 
also discussed in the physical/computational literature by [12]) enabled the rigorous iden- 
tihcation of the bell-shaped nature of the wave without [30] or with [3T] precompression. 

More recently, variants of the standard granular system in which internal oscillators or 
resonators are present in each of the lattice nodes (i.e., for each of the beads) have been 
proposed theoretically and some of them have also been realized experimentally. Arguably, 
one of the earliest examples of this type (which, however, has not yet been experimentally 
implemented, to the best of our knowledge) is the so-called cradle system [131E], where a 
local, linear oscillator was added to each bead, enabling the observation of some intriguing 
dynamical features. More recently, another variant of the granular chain, namely the 
locally resonant granular crystals or so-called mass-in-mass or mass-with-mass systems 
have been proposed and experimentally realized, respectively in [3] and m- The former 
setting involves an internal resonator within the chain (e.g. a second bead embedded in 
the principal one), while in the latter case, the resonator is external to the bead (see e.g. 
the prototype also put forth in [TH]h 

The recent experiment of [13] is expected to provide a signihcant boost to studies along 
this direction, as it offered a different so-called woodpile conhguration constructed out of 
orthogonally-stacked rods (with every second rod aligned, in this alternating 0-90 degree 
conhguration) for which it was demonstrated that the MiM/MwM description of a granular 
chain with an internal resonator can be tuned to be the relevant one [19]. In fact, such a 
system allows, depending on the parametric regime and the vibrational modes of the rods, 
to create settings where controllably two- or more internal resonator modes are relevant in 
the mathematical model description. Moreover, this experiment reported the observation 
of traveling waves for this system with persistent tails. 

The focus of the present study is to provide a rigorous analysis of the traveling waves 
in such a MiM/MwM system bearing an internal resonator. In particular, the surprising 
hnding that our analysis will rigorously establish concerns the fact that under special (i.e., 
non-generic or isolated) so-called anti-resonance conditions, the traveling waves of this 
chain will still bear a bell-shaped structure, monotonically decaying on either side. That 
is to say, for such isolated combinations of the resonator mass/wave speed, the tails of 
the wave will be absent. On the contrary, the generic setting will, in fact, involve such 
a resonant excitation through the passing the wave of the internal resonator and as such, 
it will lead to the formation of the weakly nonlocal solitary waves (nanoptera) reported 
in [13]. Our presentation of these results will be structured as follows. In section 2, we will 
provide the general setup of the problem and an integral (single component) reformulation 
thereof. In section 3, the main result will be presented and proved. Finally, in section 
4, some supporting numerical computations will be provided in order to corroborate the 
relevant results. 


MASS IN MASS SYSTEM OF GRANULAR CHAINS 


3 


2. Analysis of the problem: the model and its integral reformulation 

Following the earlier MiM/MwM studies [31HH UHl [THl [I9] , we consider the prototypical 
mathematical model of the form: 

dttX, = [(A ,_1 - X,f_, - (A, - A,+i)^] + ~k{x, - A,) 

udttXi = -k{xi - Xi) 

where p > 1. Here Aj represent the displacements of the granular beads and Xi those of 
the internal (or external) resonators, k corresponds to the (normalized) elastic constant of 
the bead-resonator coupling, while u is the (again, normalized to the mass of the beads) 
mass of the resonator particles. In the strain variables A, := Aj_i — A*, yi = Xi-i — Xi, the 
equations take the form 

,,, 3„Y, = (F4,)+ - 2{Yi)l + + k{y, - yj 

vditVi = -k[Vi - Yi) 

We hrst transform the problem to a setting which allows us to use the methods of calculus 
of variations. This is similar to the approach that we took in the earlier works [301 El] to 
treat the case of monomer chains (in the absence of resonators). 


2.1. Reduction to a single equation for $. We are now looking for traveling wave 
solutions of (jll) in the form Yi{t) = $(i — cf) and yi{t) = ^(i — cf), where we assume that <h 
will be a positive function and we set hereafter for simplicity k = 1. Plugging this ansatz 
in the system (E]) yields the following system of advance-delay differential equations 

c2$" = + 

z/c^T" = -(T - <F) 

where, we have introduced the discrete Laplacian Adiscf{x) = f{x + 1) — 2/(a;) -f f{x — 1). 

In order to restate the problem in its equivalent Fourier variable form, we introduce the 
Fourier transform and its inverse as 

/ oo poo 

fix) = / 

-OO J —oo 

The second derivative operator has the following representation 

^(0 = 


whereas 

AdiscfiO = -4sin2(7r0/(0- 

Thus, taking Fourier transform in the second equation of ([31) allows us to express 

1 


(4) 


4^(0 = 


r4>(0- 


1 — 

Plugging this last formula in the hrst equation of ([31) then yields the following equation 
for <F 

= -4sin2(7r0<h^(0 + , . , 


1 — 
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Solving for $, we obtain 

(5) $(0 


l-47r2cV^2 sin2(7r^) 

1 + U — 47r2c^Z/^2 




2.2. An integral equation for $. With the assignment A = A^, := -^/l + we can 
represent 


1 — 


= 1 


u 


1 + u — 1 + z/ — 47r^c^z^,^2 

which allows us to further rewrite (I5l) as follows 

1 


= 1 + 


( 6 ) 


*(«) = + — 


2A \27rc.^ — A 2tic^ + A/ ’ 
sin^(7r^)^ 


q2^2^2 ' 2A \27rc^ —A 27rc^ + A y 

Recall that (see [8] and also ra, EB) taking the “tent” function, A(a;) = (1 — |x|)+ or 


A(x) = 


f 1 ~ k 

X 

1 0 

X 


we have A(0 = 

Next, we compute the inverse Fourier transform of ^ . ^-kc^-a 27 rc$+A 


We have 


„27r: 


2ttc^ — A 


= 


A — x 

e c 


2nc 


sin(2;x) , 6* = "^ 

i I —^—-dz = 


z 


2c 


isgn(x). 


whence 




-1 


■^sin f^x ) sgn(a:). 


2A \27rc^ — A 27rc^ + A^ 

Thus, taking inverse Fourier transform in (|6]), we obtain the analogous to [8] (see also |12] 1 
representation for the solution of the form; 


(7) 


c2$ = A * A>P 


sin ( —X ) sgn(a:) * A * <F^. 


2Ac 


The next task is to compute the kernel 

-— sin ( —X ) sgn(a;) * A. 

2Ac \c J 

Clearly, this is an even function, being the convolution of two even functions. Thus, we 
only need to compute it for x > 0 and then we can take an even extension across zero. A 
direct computation shows 


I . [A \ , . , J ^ +-yr (sin(^)cos(^x)-sin(^x)) xe(0,1). 

2AC I 7"^ i ^ ^ 1 ^^in(^)(cos(4)-l) 


A3 


X > 1 
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At this point, we introduce a few more notations in order to rewrite ([7]) in more compact 
form. Namely, let 

^ f sin(/j) cos(/^3:)-sin(/^a:)sgn(3;) ^ ^ 

> 0, G(x) — G/j,^a(x) sin(;jx)(cos(At)-l) „g^^^^ |^| ^ ^ 

Note that the function G is compactly supported only if fi = 2n7r. 

We can now rewrite ([7]) in the form 

(8) = (A^ - 1)A * + G * <I)P. 

For positive solutions of ([H]), we can take the transformation Z = which leads us to 

(9) -1)A* Z+ G*Z = ((A^ - 1)A + G)*Z. 

Denote 

K{x) := {A^ -1)A + G{x). 

Note that K = K{A, c; x) and the problem flTT]) now reads 

(10) A^AZ^'p = Ka,c * Z. 

2.3. The anti-resonance condition. Hereafter, we restrict ourselves to the case of com¬ 
pactly supported kernel K in the convolution problem ([9]). That is, we require 

(11) /i = 2n7r, n e N 

in order to achieve that the function G (and hence K) is supported in (—1,1). We refer to 
flTTl) as the anti-resonance condition for the parameters. Notice that this is a condition 
that physically connects the mass of the resonator (or effectively the ratio of its mass to 
that of the principal bead) to the speed of the wave. We can then obtain the following 

Lemma 1. Let n G N,fi = 27rn. There exists an irrational number Aq = 1.10328..., so 
that for A > Aq, the kernel 

Ka,cA) = - l)A(a:) - 

/i 

is positive in (0,1). For A > y/2, the function K is decreasing in (0,1). 


Proof. Since K is obviously an even function, it suffices to consider the case 0 < x < 1. 
Taking the derivative of K, we see that 

K'{x) = — cos(/ix) — (A^ — 1) < 1 — (A^ — 1) < 2 — A^ < 0 


so long as A > y/2. Thus, K is decreasing in (0,1) for A > \/2. 

Next, we study the positivity of K. Again, it suffices to consider x G (0,1). We have 


N / a 2 \ sin(ux) _ , 

K{x) = (A^ — 1)(1 — x) -= (1 — x) 

T 


A' 


sin(/ix) 
/i(l -x) 


Since = 2n7r, we can rewrite the last expression as follows 


K{x) = (1 — x) 


A^ - 1 + 


sin(27m(l — x)) 
27rn(l — x) 
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Now, since 

^ sin z ^ sin z 

—oo = inf -= inf -= —0.217234...., 

^G(0,27r) Z 2G(0 ,oo) Z 

we conclnde that K{x) > 0, provided A> Aq = yjl + oq = 1.10328... 


□ 


We can now proceed to establish onr main resnlt. 


3. Main result: Bell-shaped Traveling Waves persist Under Suitable 

Anti-Resonance Conditions 


We have the following main resnlt. 

Theorem 1. Let ^ = 2mi for some integer n = 1, 2,.... 

Then, for A > \/2, the equation (1T0|) (and hence (jS]) j has bell-shaped solution T. That 
is, there A : R ^ R+ a positive, even, C°° smooth function, so that $ is decaying in 
(0, oo). In addition, $ is doubly exponentially decaying, just as the travelling waves in the 
classical monomer case. 

For \/2 > A > Ao = 1.10328..., the problem (and hence ([5]) ) has a positive solution 
<h, which we cannot guarantee to be bell-shaped\j. 

We perform the proof in several steps, nnder the different assnmptions on A. Following 
the idea in [301 lOT] , we let g = 1 ^ and we set np the following maximization problem 

. . / {K*z,z) ^ max 

\ /j^i |2;(a:)|''dx = 1 

Note that q G (1,2). The plan of actions is as follows: we hrst show that (IT2|) has a 
solntion with the desired properties. After that, we derive its Enler-Lagrange eqnation. 
This is of conrse closely related to (fTOj) . except for the Lagrange mnltipliers, which need 
to be adjusted. 


3.1. Existence for the maximizers. We will show that the following lemma holds. 

Lemma 2. Let A > Aq. Then, the maximization problem flTO]) has a solution. Moreover, 
this solution is positive. If in addition A > \/2, then it is also bell-shaped. 


First, we show that the quantity {K * z, z) is bounded, if z satishes the constraint. 
Indeed, by Holder’s and Young’s inequalitjH 

\{K*z,z)\ < \\K*z\\l,'\\z\\l, < \\K\\^^\\z\\l„ 

q 

But K is a bounded function with support in (—1,1), hence K G L 29-2 . Next, we show 
that an eventual solution to flT^ is necessarily positive. Indeed, for any function z, we 
have that the function w := \z\ satisfies the constraint and moreover 


(K * z, z) 


K{x — y)z{x)z{y)dxdy < 


K{x — y)w{x)w{y)dxdy = {K *w,w). 


^The latter feature will, however, be illustrated in our case example numerical computations of the next 
section. 

^Note that since g = 1 -|- i < 2, we have that .^q-i ^ 
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Thus, the supremum in ffT^ may be taken over all 2 ; > 0. In fact, one can see that if 2 ; is 
not non-negative, then w = \z\ provides a bigger value and hence 2 ; may not be a solution 
to (|T^ . Denote 

jmax ^ if Z^Z). 

lhllL'2=l 

Clearly, > Q. Pick a maximizing sequence, say Zn'- Zn> 0. That is H^nllL-? = 1 and 


{K*Zn,Zn)^J'^^^. 

We will show that (an appropriate translate of) Zn converges strongly (in L^) to a solution 
Zq. To that end, we apply the concentration compactness method of Lions, [23]. The 
outcome is that (after we pick a subsequence, which we call again Zn), one of three scenarios 
occur: 

• (tightness) There exists G R^, so that for each e > 0, there is i? = R{e), so that 

ryk+R 

/ z\{x)dx > 1 — e 
Vk-R 


(vanishing) For every R> 0, there is 

py^R 


lim sup 


k^oo 


zl{x)dx = 0. 


y Jy-R 


(dichotomy) There exists a G (0,1), so that for any e > 0, there are R = R{e) and 
Rk — 00 , i/fc G so that for all large enough fc. 


rVk+R 


r 

/ z\{x)dx — a 

A 

rh 

•0 

/ z\{x)dx — (1 — ck) 

Jyk-R 


J\y-yk\>Rk 


< e'^ 


We will proceed to show that vanishing and dichotomy may not occur, which will leave 
us with tightness. In the tightness scenario, we will easily show that a translate of Zk will 
converge strongly to Zq. 

Vanishing does not occur: 

Assume that it does. Let e > 0, R = 10 and ko = ko{e) is so large that 


(13) 


fy+10 


sup 


zl{x)dx < e^. 


y Jy-W 

for all k > ko. We have for each |/o G R^, 

ryo+5 

/ K * Zk{y)zk{y)dy < ||zfc||x,<j(j^Q_5^j^Q+5)|| JF * ^fc||x,9'(yp_5^yp+5)- 

dyo-5 

Since suppK C (—1,1), and the integration is over (yo — 5,yo + 5), it follows that 

X(yo—5,yo+5)R * ^k{y') — X(j/o—5,j/o+5)A" * [ZkX{yo—6,yo+G)]- 

We have by Young’s inequality that 
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It follows that 

/•yo+5 

/ K* Zk{y)zk{y)dy < C\\zk\\l,^y^_io^y^+^o) < Ce^-<^\\zk\\l,^y^_^Q y^^^Qy 

dyo -5 

where in the last inequality, we have used fflSj) . Summing the last inequality over y^ = 
0 , ±1 ,... yields 

/ yo +5 

K * Zk{y)zk{y)dy < Cie^~'‘\\zk\\L<i = Cie^~‘>, 

for all large enough k and for all e > 0. On the other hand {K * Zk, Zk) which will 

be a contradiction, if we have had selected e small enough. Note that again, we have used 
q < 2. Thus, vanishing may not occur. 

Dichotomy does not occur: 

Assume that it does. Let e > 0 and denote 


4iy) ■= Zk{y)X(y,-R,y,+R){y); 44) = ^kiy)X\y-y,\>R4y)^ 

so that \\zk — z\ — = 0(e). From support considerations and the Young’s estimates 

that we have provided earlier, it is clear that for all large enough k 

{K * Zk, Zk) = {K* zl zl) + {K* zl zl) + 0(e). 


Now, from the dehnition of 

(^*4,4) = 114114^ * 

Similarly, 



7^ 

\ ^ .-^1 l|2 

Tj-iTj-; A J WZkWq 
W^kWQ 


J“""(a2 + 0(e9)). 


(iF*zL4)<^““((l 

Putting everything together yields 

{K*Zk,Zk)<J^^^4a^ + {l- 


-af + 0{e^)). 
a)2 + 0(e'?)) + 0(e). 


This is again produces a contradiction (with a judiciously small choice of e), since 
{K * Zk, Zk) —)■ Q!^ + (1 — aY < 1. Thus, dichotomy does not occur either. 

Thus, tightness is the only alternative. Now, consider Zniy) ■= Zniy — y-n)- Note that Zn 
satishes the constraint and also {K * Zn, Zn) = {K * Zn, Zn) We have that for all 

e > 0, there exists R = R{e), so that 

(14) [ zl{y)dy > 1 - e. 

J-R 

for all n. Now, by weak compactness in L^, it follows that Zn has a weakly convergent 
subsequence (denoted again by Zn for convenience), with a limit say zq. 

We now claim that the sequence {K * Zn}n is strongly pre-compact in . Indeed, by 
the Young’s inequality, we have for all n 

ll-Y * Zn\\]Yl,q' < ||/F * -|- \\K * ZnWiq' — ^IIIknllL-? = C*|| A'H, 
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1 9 

where by inspection K eW ’ 29 - 2 . Also 

[ \K* i„{y)\-/dy < ||;r||«; 

J\y\>R+l 

By the compactness criteria in L'^' spaces (i.e. the Riesz-Tamarkin condition), it follows 
that {K * Zn}n is pre-compact. This, together with the pointwise limit K * Zn ^ K * zq, 
which follows from the weak convergence of Zk implies that for some subsequence 


We now have by Holder’s 


K *Zni- K* Zolli,' ^ 0. 


\{K * Zni,Znt) - {K * Zo,Zo)\ < 

< 


\{K * Zni - K * Zo, Zni)\ + \ (^ 0 , K * Z^ - K * Zo)\ 
* Zni — K * ^oIIlT (II^oIIl'J + ll-Snillis) —)■ 0. 


Thus, 


{K*zo,zo) = J^^^. 


On the other hand, by the lower semicontinuity of the norm with respect to the weak 
convergence ||. 2 o||l'? < 1- But then ||. 2 o||l 9 = 1; since otherwise 


J-^-^ = {K*zo,Zo) = \\zo\\UK* 


Zo 


Zo 


IkolU.’lkoll 


L<i 


-) < J^nUoWh < J 


a contradiction. Thus, we have shown that the limit zo is indeed a solution to the maxi¬ 
mization problem (I12p . 


3.2. Bell-shapedness of the solution in the case A > \/2. In this case, we have that 
the kernel K is bell-shaped. We show now that the solution zq is in addition bell-shaped. 
In order to explain the setup, we need a few dehnitions. 

For a function / ; R^, one dehnes the distribution function 


df{a) := meas{x : |/(a:)| > a}. 

Clearly, the function df{a) is non-increasing, whence one can dehne its “inverse” as follows 

f*{t) = inf{s : df{s) < t}. 

We call /* ; Ri'|_ —?■ R^ the non-increasing rearrangement oi f. Note that the two functions 
have the same distribution function, that is df{a) = df*{a), so in particular ||/||i,p = 
11/11 LP : 0 < p < oo. Finally, dehne the even function f*{t) = /*(2|t|), which also satishes 
II/IIlp = 11/^ I Up : 0 < p < oo. Clearly, a function is bell-shaped if and only if /^ = /. In 
this setting, we have the Riesz rearrangement inequalitjH, [1], , which states 



y)g{y)h{x)dxdy < 



y)g*{y)h*{x)dxdy. 


Using that fact and since K is bell-shaped, it is clear that in the maximization problem 
m, we can restrict our attention to bell-shaped entries .s. Indeed, taking an arbitrary 


3 


The original inequality appeared implicitly in the work of Riesz. 
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function z, note that would satisfy the constraint = \\z\\Lq = 1 and moreover, 

by the Riesz rearrangement inequality 


{K*z,z)= / K{x — y)z{y)z{x)dydx < / K{x — y)z*{y)z*{x)dxdy = {K * z*, z*). 
Jn? Jn? 

Based on the above formulation, one may derive the existence of Zq the same way as before. 
Note however that one may completely circumvent the Lions concentration compactness 
arguments in the bell-shaped case. Indeed, assuming that 2 ; is bell-shaped and satishes the 
constraint \\z\\Lq = 1, we have for every x > 0, 

f*X f*00 

\z{x)\^x< / \z{y)\^dy< / \z{y)\^dy = 1/2. 

Jo Jo 

Hence \z{x)\ < 2~^x~^l^. Thus, any maximizing sequence {zn} of bell-shaped functions will 
prodnce a pre-compact in L'^ sequence K*Zn- Indeed, we have, similarly to the argnments 
above that K * Zn & . In addition, by support considerations for any R > 10, we have 

for every n 

f \K * Zn{x)\‘^ dx < C f \z{x)\‘^ dx < C f . dx < CR^~~ 0, 

J\x\>R J\x\>R-l J\x\>R-l\^r '^ 

as i? —)■ 00 , since q < 2 < q'. Thus, the elementary pointwise bounds obtained from 
bell-shapedness help us verify immediately the Riesz-Tamarkin criteria and hence the pre¬ 
compactness oi K * Zn in , without having to resort to the Lions’ theory. After that, we 
finish by standard arguments as before. 


3.3. The Euler-Lagrange equation for flT2|) . The next step is to derive the Euler- 
Lagrange equation that the solution zq of flT^ satishes. We proceed as follows. Take 
e G and a test function h. Since Zq is a maximizer for flT2]) . it must be that 


(K* 


Zq + eh 

Iko + ^h\\lq ’ 


zp + eh 
Iko + ^hWiq 


for all e, h. We can rewrite this as follows 


( 15 ) {K * (zq + eh), Zq + eh) < J'^^^Wzq + ^h\\\q 

Taking Taylor expansions in e, we hnd 

{K * {zq + eh),ZQ + eh) = {K * zq, zq) + 2e{K * ZQ,h) + 0{e^) = 

= r^^^ + 2e{K*ZQ,h)+0{e^), 

\\zo + eh\\\q = {1 R eq j zl~^{x)h{x)dx + 0{e^)Y = I+ 2eq{zl~^,h) + 0{e^). 

Plugging in this in flT^ . we obtain 

2{K * zo, h) = 2J^^^q{zl-\ h) + 0(e), 
which is valid for all e, h. Taking limit as e —?■ 0, we obtain 

{K * Zq- h) = 0 - 
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Since this is satisfied for all test functions h and g — 1 = we obtain K * zq = qz^. If we 
now take zo = fiZ, 

K * Z = q^p~^Zp. 

Clearly, with an appropriate choice of /i, we can arrange so that Z indeed satisfies (ITU]) . 

We now proceed to show that Z and <h = are both C°°. First, let us see that from 
the equation flT 0 |) . the function Z never vanishes. Indeed, assume the opposite, Z(xo) = 0. 
Then, 

0 = A^c^Z^/P{xo) = J K{y)Z{x - y)dy. 

Since K{y) > 0 in (—1,1), it would follow that Z\(^xq-i,xo+i) = 0. Iterating this argu¬ 
ment yields that Z{x) = 0 for all x, a contradiction. Now, we have from (fTOj) that <h is 
differentiable and moreover 

= K' * 

It follows that is continuous. By induction, we can show 

^ K' 

and hence is continuous, since < I are all continuous and <F 7 ^ 0. For integer 

values of p, this is obvious. Even for non-integer values of p, is as smooth as <F, since 
<h(x) > 0. 

Finally, let us establish the double exponential rate of decay for the bell-shaped solutions 
of (|5|). So, assuming A > \/2, we have that $ is bell-shaped and K is decaying for x > 0. 
We have 

rx+l 

A^c^<h(x) = / K{y)^^{x — y)dy. 

Jx-l 

Denote l> ;= ^ ^ so that 

= W ^\— / K{y)^P{x - y)dy. 

|| l°c Ja:-l 

Using the bell-shapedness of both K and <h we conclude 

l>(x) < $^(x — 1 ). 

Iterating this inequality for n < x, 

<l>(x) < |)^"(x — n). 

which yields the double exponential rate of decay for $ and hence for $. This completes 
the proof of Theorem [TJ 

In order to corroborate the above conclusions, as well as to examine the cases which do 
not adhere to the conditions of the above theorem (both cases of anti-resonance, but with 
A< V2, as well as cases where the anti-resonance condition is not satisfied), we turn to 
some case example numerical computations in the next section. 
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4. Numerical results 


4.1. Anti-resonance scenario (/i = 2?™). When jj, = 27rn, n G N, we obtained the 
expression of kernel (as shown in Lemma 1) in the convolution equation: 

= Ka,c * (<h)P. 


When a solution of <h is obtained, one can similarly hnd 4/ from: 


(16) 


^(0 


_1_ sin^«) 

1 + U — 




and its corresponding integral equation 


(17) ^ = Ma,c * 

where 

MaA^) = ~ ^ —^sgn(2:)]X(-i,i)(a;). 

Using the above expressions, the solution for both the bead and the resonator strains can 
be not only theoretically analyzed but also numerically computed. 

Due to the similarity between Ma^c and /’7 a, c, we will mainly discuss the properties 
of Ka,c and those of Ma,c directly follow. Since Ka,c is an even function and has hnite 
support (—1,1), we only need to consider its properties on (0,1). Noting that = 

= —A^, /7(^^(0“) = we appreciated that /7 a, c is continuous 

everywhere but not continuously differentiable at x = —1,0,1. We now separate the 
different cases of interest. 

(1) The case analyzed fully in the previous section (i.e., the “best case scenario”): 
when A > y/2, K'j^J^x) < 0 for x G (0,1) so that Ka,c is increasing on [—1,0] and 
decreasing on [0,1]. This regime also implies the positivity of /7 a,c since /7 a,c(—1) = 
//a,c( 1) = 0. Due to these good features of KaAAi and as a by-product of the 
proof of the previous section, our numerical computations conhrm that the system 
bears bell-shaped solutions for $ and T, as illustrated in the typical results shown 
in Fig. [H 

(2) The “intermediate” case: when 1.10328 Aq < A < \/2, Ka,cA) ^ 0 still holds 
for all X but the kernel is no longer monotone on (—1,0) or (0,1). Although in 
this category KaAA has non-monotonic variations (with the number of local min¬ 
ima being n), the corresponding solution of <F still features a bell-shaped prohle 
as it did in the previous case (see Fig. [2] for such examples in the case of differ¬ 
ent anti-resonances). The difference between the (non-mo notonic) kernel and the 
(monotonic) solution is, apparently, caused by the smoothing effect of convolution. 

(3) Turning now to cases for which our analytical arguments can no longer be used 
(even in a qualitative way), we hrst consider: 1 < Ai^„ < A < Aq. Here, KaAA is 
neither necessarily positive, nor decreasing on (0,1). However, as Fig. E] suggests, 
we notice from the numerical results that a bell-shaped solution of <h and 4/ could 
still exist as long as A is not too small, i.e. A > Ai^„ > 1. Here Ai^„ decreases over 
n and it is the threshold above which our numerical method is able to converge to 
a solution, as shown in the right panel of Fig. 01 
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Figure 1. The top panels show the kernel of the convolntion equation 
Ka,c{x) and corresponding solution $ and T with A = 1.5 and /i = 27r. 
In the middle panels and bottom panels, /x is set as dvr and Gtt, respectively, 
corresponding to higher anti-resonances. 
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X X 




Figure 2. This figure is structured similarly as Fig. [T]but A has changed 
to 1.2 instead of 1.5 here. The top panels show the kernel of the convolution 
equation Ka,c{,x) and corresponding solution $ and T with A = 1.2 and 
/i = 27r. In the middle panels and bottom panels, /x is set as dvr and bvr, 
respectively. 
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Figure 3. This figure is similar to the above figures and the only difference 
is that A = 1.09 here. The top panels show the kernel of the convolution 
equation Ka,c{x) and corresponding solution $ and T with A = 1.09 and 
/i = 27r. In the middle panels and bottom panels, p is set as dvr and Ovr, 
respectively. 
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(4) When 1 < A < Ai^^i the substantial intervals of negative values of Ka,c{x) lead 
our numerical scheme to failure of convergence to any solution for $ and 4/. The 
plot of a typical example of Ka,c{,x) within this category is given in the left panel 
of Fig. m The dependence of Ai^n on n is shown in the right panel of the hgure, as 
indicated above. 




Figure 4. The left panels shows the kernel of the convolution equation 
Ka,c{x) with A = 1.02 and fi = 2 ti while how Ai^n decreases over n is 
revealed in the right panel. 

(5) Resonance case with /i ^ 27 rn . In this case, the conditions for calculating the 
Fourier transform on $ and 4/ fail. By utilizing Fourier series instead (on a hnite 
computational interval), we obtain that Ka,c = — 1) max (1 — |a;|, 0) + G{x) has 

non-decaying oscillatory tails on both wings. Moreover, despite the fact that Ka,c 
is neither increasing nor non-negative on (—cxo, 0), we are able to obtain numerical 
solutions for <|) and T from either the convolution equations (i.e. the integral 
equations) or the system of advance-delay differential equations, as shown in Fig. |5l 
It is interesting that the solution also has a bell shape in the center, but it possesses 
oscillatory tails on both wings. The computations of this resonant case and their 
subtleties (from a numerical perspective), as well as the different methods utilized 
to obtain the solutions are analyzed in detail elsewhere [53]. The rigorous analysis 
which is the main emphasis of the present work cannot, unfortunately, presently 
provide dehnitive insights about the latter case. This remains an important open 
problem for future study. 

5. Conclusions & Future Challenges 

In the present work, we have explored systems in the form of Mass in Mass or Mass 
with Mass dynamical lattices possessing an internal resonator, studying their traveling 
waves. In particular, we could explore the case of so-called anti-resonances whereby for 
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Figure 5. The left panel shows the kernel of the convolution equation 
Ka,c{x) and the right shows the corresponding solution $ and T with A = 1.2 
and fi = 2% — 0.5. 

a particular set of “quantized” relations between the resonator mass and the wave speed 
(and for suitably small resonator masses), a bell-shaped traveling wave could be rigor¬ 
ously proven to exist. Interestingly, the bell-shaped waves were numerically found to exist 
when the anti-resonance condition applies, even beyond the mass threshold for which our 
proof holds. However, a secondary threshold was identihed for sufficiently large resonator 
masses, beyond which our iterative scheme for the numerical identification of the waves no 
longer converged. Finally, when the anti-resonance condition is not upheld, typical hndings 
suggest the existence of a wave of non-vanishing tail, in line also with recent experimental 
observations [18] . 

Naturally, many directions of future research emerge from the present (and recent) con¬ 
siderations of this novel class of systems. From a rigorous perspective, understanding the 
phenomenon of resonances and the formation of traveling waves with tails would be es¬ 
pecially interesting. Equally interesting from an experimental perspective appears to be 
the actual experimental setting whereby an initial excitation of the chain, by construction, 
produces such tails which are arising only on one wing but not the other. Additionally, as 
indicated in [18] , a distinct experimental possibility is that of bearing multiple resonators 
rather than one in the context of the woodpile configuration. Understanding how one 
vs. more such resonators may affect the observed phenomenology is another important 
direction for both theoretical and experimental future studies. 
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